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Abstract 



o 

. For the pp and pp systems, we have used all of the extensive data of the Particle 

\ Data Group[K. Hagiwara et al. (Particle Data Group), Phys. Rev. D 66, 010001 

(2002).]. We then subject these data to a screening process, the "Sieve" algorithm [M. 
M. Block, physics/0506010,], in order to eliminate "outliers" that can skew a fit- 
r~| , With the "Sieve" algorithm, a robust fit using a Lorentzian distribution is first made to 

all of the data to sieve out abnormally high Ax?, the individual i*^ point's contribution 
to the total x^- The x^ fits are then made to the sieved data. We demonstrate that 



we cleanly discriminate between asymptotic In s and In'^ s behavior of total hadronic 



X 

I cross sections when we require that these amplitudes also describe, on average, low 

energy data dominated by resonances. We simultaneously fit real analytic amplitudes 
to the "sieved" high energy measurements of pp and pp total cross sections and p- 
values for > 6 GeV, while requiring that their asymptotic fits smoothly join the 
the app and cr^p total cross sections at y/s =4.0 GeV — again both in magnitude and 
slope. Our results strongly favor a high energy In'^ s fit, basically excluding a Ins fit. 
Finally, we make a screened Glauber fit for the p-air cross section, using as input our 
precisely-determined pp cross sections at cosmic ray energies. 



1 Introduction 



This paper consists of 3 parts: 

1. The introduction of the "Sieve" algorithm of M. M. BlockpQ, used to "sieve" out 
"outhers" from the Particle Data Group [TU] compilations of cross sections and p— values 
(the ratio of the real to the imaginary portion of the forward scattering amplitude). 

2. A new, constrained method |7j of fitting p and pp cross sections and p— values, showing 
that the Froissart bound is saturated. This work was done in conjunction with Francis 
Halzen. 

3. A Glauber calculation for the p-air production cross section, using screening, an un- 
published work done in collaboration with Ralph Engel. 

2 Part 1: The "Sieve" Algorithm 

We now outline the adaptive Sieve algorithm pP that minimizes the effect that "outliers" — 
points with abnormally high contributions to — have on a fit when they contaminate a data 
sample that is otherwise Gaussianly distributed. Our fitting procedure consists of several 
steps: 

1. Make a robust fit of all of the data (presumed outliers and all) by minimizing Aq, the 
Lorentzian squared with respect to a, where 

N 

Alia; x)=Y.^n{l + 0.179Ax^(x,; a)} , (1) 

i=l 

with a = {«!,..., ajv/} being the M-dimensional parameter space of the fit. x = 
{xi, . . . , xn} represents the abscissa of the experimental measurements y = {yi, . . . , i/n} 
that are being fit and Axiixf, a.) = ^Mizy^ii^^ jg individual contribution of 

the z*^ point, where y{xi; a) is the theoretical value at Xi and cTj is the experimental 
error. It is shown in ref. P that for Gaussianly distributed data, minimizing Aq gives, 
on average, the same total Xmin = YliLi ^xK^u ^) from eq. ((T)) as that found in a 
conventional fit, as well as rms widths (errors) for the parameters that are almost 
the same as those found in a fit- 

A quantitative measure of whether point i is an outlier, i.e., whether it is "far away" 
from the true signal, is the magnitude of its Axi{xi;ot) = (^Vizy^^i^^!!^ _ xhe reason 
for minimizing the Lorentzian squared is that this procedure gives the outliers much 
less weight w in the fit {w oc 1/ y^Axf (x^; a) ), for large Axf (x^; a)) than does a x^ 

fit {w oc ^ Axiixi] a) ), thus making the fitted parameters insensitive to outliers and 
hence robust. For details, see ref. [1 . 

If Xmin is satisfactory, make a conventional x^ fit to get the errors and you are finished. 
If Xmin is ^ot satisfactory, proceed to step El 
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2. Using the above robust Aq fit as the initial estimator for the theoretical curve, evaluate 
Axfixf, a), for each of the experimental points. 

3. A largest cut, Axfixf, a) max, must now be selected. We start the process with Axiixf, a) 
9. If any of the points have Axfixi] a) > Axfixi] Q:)jnax, reject them — they fell through 
the "Sieve". The choice of Axfixf, a)inax is an attempt to pick the largest "Sieve" size 
(largest Axfixi] Q:)max) that rejects all of the outliers, while minimizing the number of 
signal points rejected. 

4. Next, make a conventional fit to the sifted set — these data points are the ones that 
have been retained in the "Sieve". This fit is used to estimate Xmin- Since the data 
set has been truncated by eliminating the points with Axfixf, a.) > Axfixi] Q:)max, we 
must slightly renormalize the Xmin found to take this into account, by the factor TZ. 
For Axj^i^^ = 9, 6, and 4, the factor 71 is given by 1.027, 1.140 and 1.291, whereas the 
fraction of the points that should survive this cut — for a Gaussian distribution — is 
0.9973, 0.9857 and 0.9545, respectively. A plot of 7l~^ as a function of ^ax given 
in Figure ^ which is taken from ref. jT]. 

If the renormalized Xmin^ ^-^v x Xmin is acceptable — in the conventional sense, using 
the ordinary distribution probability function — we consider the fit of the data to the 
model to be satisfactory and proceed to the next step. If the renormalized Xmin is not 
acceptable and Axfixf, Q:)max is not too small, we pick a smaller Axjixf, Q;)max and go 
back to stepEl The smallest value of Axfixf, Q:)max that we used is Ax^ixf, Q;)max = 4. 

5. From the fit that was made to the "sifted" data in the preceding step, evaluate the 
parameters a. Next, evaluate the M x M covariance (squared error) matrix of the 
parameter space which was found in the x^ fit- We find the new squared error matrix 
for the A^ fit by multiplying the covariance matrix by the square of the factor r^2. 
From Figure HJ we find that r^2 ~ 1.02,1.05 and 1.11 for Axfixf, a.)uia.x = 9, 6 and 
4, respectively . The values of r^2 > 1 reflect the fact that a x^ fit to the truncated 
Gaussian distribution that we obtain — after first making a robust fit — has a rms (root 
mean square) width which is somewhat greater than the rms width of the x^ fit to the 
same untruncated distribution [Tj. 

The application of a x^ fit to the sifted set gives stable estimates of the model parameters 
oc, as well as a goodness-of-fit of the data to the model when Xmin is renormalized for the 
effect of truncation due to the cut Axfixi] Q:)max- One can now use conventional probabilities 
for x^ fits, i.e., the probability that x^ is greater than 71 x Xmin^ fo^ the number of degrees of 
freedom u. Model parameter errors are found by multiplying the covariance (squared error) 
matrix of the conventional x^ fit by the appropriate factor [r^i)"^ for the cut Axfixf, Q:)max- 

3 Part 2: Saturating the Froissart Bound 

High energy cross sections for the scattering of hadrons should be bounded by a ~ In^ s, 
where s is the square of the cms energy. This fundamental result is derived from unitarity 
and analyticity by Froissart [2], who states: "At forward or backward angles, the modulus 
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of the amplitude behaves at most hke s In^ s, as s goes to infinity. We can use the optical 
theorem to derive that the total cross sections behave at most like In^ s, as s goes to infinity" . 
In this context, saturating the Froissart bound refers to an energy dependence of the total 
cross section rising no more rapidly than In^ s. 

The question as to whether any of the present day high energy data for pp and pp cross 
sections saturate the Froissart bound has not been settled; one can not unambiguously 
discriminate between asymptotic fits of Ins and In^ s using high energy data onlyOHj. We 
here point out that this ambiguity is resolved by requiring that the fits to the high energy data 
smoothly join the cross section and energy dependence obtained by averaging the resonances 
at low energy. Imposing this dualityj^j condition, we show that only fits to the high energy 
data behaving as In^ s that smoothly join (in both magnitude and first derivative) to the low 
energy data at the "transition energy" (defined as the energy region just after the resonance 
regions end) can adequately describe the highest energy points. This technique has recently 
been successfully used by Block and Halzen|6^ to show that the Froissart bound is saturated 
for the 7p system. 

We will use real analytic amplitudes to describe the data. The total cross sections a 
are found from the optical theorem and p is the ratio of the real to the imaginary portion 
of the forward scattering amplitude. As shown in ref. [7j, in the high energy limit where 
s 2mu, we can write cr^ and p^, along with the cross section derivatives ^^J^^^ , as sums 
and differences of even and odd amplitudes, i.e., 

.± =c„ + c.lnm+c,l„^f^)+,3,,(^)""± ^ii)"". (2) 

m/ Vm/ Vm/ Vm/ 

f^ = ^{^c, + c.m (^) - p,, cot f ^) (^y' + l^/,(o) 



Cl 




±5tan — - , (3) 



a[[v/m) y[v/m)j [ [u/mj J J 

±6{{a-mu/m)r''}, (4) 

where the upper sign is for pp and the lower sign is for pp scattering. The exponents fi and a 
are real. The real constant /+(0), appearing only in the p- value, is the subtraction constant 
at = needed to be introduced into a singly-subtracted dispersion relation [Hj, [HI- We note 
that eq. is linear in the real coefficients cq, ci, C2, and 6, convenient for a fit to the 
experimental total cross sections and p-values. Throughout we will use units of u and m in 
GeV and cross section in mb, where m is the proton mass. 
It is convenient to define, at the transition energy uq, 

2 

Co + Ciln(z/o/m) + C2 ln^(z/o/m) + P-p, {uq / mY'^ , (5) 
(7+(z/o/m) - a~{uo/m) 
2 

5{uo/mr-\ (6) 



Aa 
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1 / da+ da- 



2 \d{i'/m) (i(z//m) j 



Cl 



1 1 [2 ln(z/o/m) 



Am 



(icT^ da" \ 



u=uo 



2 \d{u/m) d{i'/m) 
= 6{{a-l){uo/mr-'}. {8] 
Using the definitions of (Xav, Act, and Am, we now write tlie four constraint equations 



(z/o/m) 



2 ln(z/o/m) 



"iav - Cl < -— > - C2 



(9) 



Co = o-av-ciln(z/o/m) -C2 In (z/o/m) -/3p/(z/o/m)^ , (10) 

« = l + ^(z/o/m), (11) 

S = A(T(z/o/m)i~", (12) 

that utihze the two slopes and the two intercepts at the transition energy z/q, where we join 
on to the asymptotic fit. We pick uq as the (very low) energy just after which resonance 
behavior finishes. We use /i = 0.5 throughout, which is appropriate for a Regge-descending 
trajectory. In the above, m = rrip is the proton mass. 

Our strategy is to use the rich amount of low energy data to constrain our high energy fit. 
At the transition energy z/q = 7.59 GeV, corresponding to a cms (center of mass) energy of 
a/sq = 4 GeV, the cross sections cr+ (z/o/m) and (T~(z/o/m), along with the slopes (^^(^7^) 

and ( d{{u/m) ) ' ^'^^ used to constrain the asymptotic high energy fit so that it matches 
the low energy data at vq. We picked z/o much below the energy at which we start our high 
energy fit, but at an energy safely above the resonance regions. Very local fits are made 
to the region about the energy in order to evaluate the two cross sections and their two 
derivatives at that are needed in the above constraint equations. We next impose the 
4 constraint equations. Equations Q, (fTUj) . (fTTj) and (fT^ . which we use in our fit to 
Equations El and El For safety, we start the data fitting at an energy z/min = 18.25 GeV, 
corresponding to the cms energy, 1/5^^;^ = 6.0 GeV, appreciably higher than the transition 
energy. 

We stress that the odd amplitude parameters a and 5 and hence the odd amplitude 
itself, 5(z//m)"~^, are completely determined hj the experimental values Am and Aa at the 
transition energy uq. Thus, at all energies, the differences of the cross sections a" — cr+ (from 
the optical theorem, the differences in the imaginary portion of the scattering amplitude) 
and the differences of the real portion of the scattering amplitude are completely fixed before 
we make our fit. Further, for a In^ s (Ins) fit, the even amplitude parameters Cq and Pp 
are determined by Ci and C2 (ci only) along with the experimental values of cTav and mav at 
the transition energy uq. In particular, for a In^ s (Ins) fit, we only fit the 3 (2) parameters 
Cl, C2, and /(O) (ci and /+(0)). Since the subtraction constant /+(0) only enters into the 
p-value determinations, only the 2 parameters Ci and C2 of the original 7 are required for a 
In^ s fit to the cross sections cr^, which gives us exceedingly little freedom in this fit — it is 
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indeed very tightly constrained, with not much latitude for adjustment. The cross sections 
cr^ for the Ins fit are even more tightly constrained, with only one adjustable parameter, Ci. 

Table^summarizes the results of our simultaneous fits — using the 4 constraint equation — 
to the available accelerator data from the Particle Data Group |^ for cipp, app, ppp and ppp, 
after using the "Sieve" algorithm. Two ^xlmsix cuts, 6 and 9, were made for ln^(z//m) fits. 
The probability of the fit for the cut Axf^;^^ = 6 was ~ 0.2, a very satisfactory probability 
for this many degrees of freedom, and we chose this data set. As seen in Tabled the fitted 
parameters are very insensitive to this choice. The same data set {^Ximax ~ ^ ^^t) 

Table 1: The fitted results for a 3-parameter fit with a ~ In^ (z^/mp) and a 2-parameter fit with 
a ~ In^v/ nip) to the total cross sections and p- values for pp and pp scattering. The renormahzed 
X^/'^minj taking into account the effects of the ^Xiraax ^^^^ given in the row labeled TZ x Xmin/^- 
The errors in the fitted parameters have been multiplied by the appropriate r^2- The proton mass 
is nip and the laboratory nucleon energy is ly. 





a ~ In^ (z^/mp) 


a ^ ln(i'/mp) 


Parameters 


■^i max 


-^t max 




6 


9 


6 


Even Amplitude 


Co (mb) 


37.32 


37.25 


28.26 


ci (mb) 


-1.440 ±0.070 


-1.416 ±0.066 


2.651 ±0.0070 


C2 (mb) 


0.2817 ±0.0064 


0.2792 ± 0.0059 




/3p/ (mb) 


37.10 


37.17 


47.98 




0.5 


0.5 


0.5 


/(O) (mb GeV) 


-0.075 ±0.59 


-0.069 ±0.57 


4.28 ± 0.59 


Odd Amplitude 


6 (mb) 


-28.56 


-28.56 


-28.56 


a 


0.415 


0.415 


0.415 


Amin 


181.6 


216.6 


2355.7 


^ ^ Xmin 


201.5 


222.5 


2613.7 


U (d.f). 


184 


189 


185 


^ X xLJ^ 


1.095 


1.178 


14.13 



also used for the ln(z//m) fit. The probability of the ln(i^/m) fit is << 10^^^ and is clearly 
ruled out. This is illustrated graphically in Fig.|2K- 

We note that when using a ln^(z//m) fit before imposing the "Sieve" algorithm, a value 
of x^/d.f.=5.657 for 209 degrees of freedom was found, compared to x^/d.f.=1.095 for 184 
degrees of freedom when using the Ay? , = 6 cut. The "Sieve" algorithm eliminated 25 
points with energies ^/s > 6 GeV (5 app, 5 (Xpp, 15 Ppp), while changing the total renormahzed 

from 1182.3 to 201.4. These 25 points that were screened out had a contribution of 
980.9, an average value of 39.2. For a Gaussian distribution, about 3 points with Ax^ > 6 
are expected, with a total contribution of slightly more than 18 and not 980.9. We see 
that the "Sieve" algorithm has rid the two data sets of outliers. 
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Figure Et^) shows the individual fitted cross sections (in mb) for pp and pp for ln^(z//m) 
and ln(i//m) for the cut ^xlraa.^ = 6 in Table HJ plotted against the cms energy, i/i, in GeV, 
and Figure |2b) shows the individual fitted p-values for pp and pp plotted against the cms 
energy, -y/i. The data shown are the sieved data with i/i > 6 GeV. The ln^(z//m) fits to 
the data sample with Ax^^^^^ = 6, corresponding to the solid curve for pp and the dash- 
dotted curve for pp, are excellent, yielding a total renormalized = 201.5, for 184 degrees 
of freedom, corresponding to a fit probability of ~ 0.2. On the other hand, the ln(z//m) 
fits to the same data sample — the long dashed curve for pp and the short dashed curve for 
pp — are very bad fits, yielding a total = 2613.7 for 185 degrees of freedom, corresponding 
to a fit probability of << 10~^^. In essence, the ln(z//m) fit clearly undershoots all of the 
high energy cross sections. The ability of nucleon-nucleon scattering to distinguish cleanly 
between an energy dependence of ln^(z//mp) and an energy dependence of ln(z//mp) is quite 
dramatic. 

A few remarks on our ln^(z//m) asymptotic energy analysis for pp and pp are in order. 
It should be stressed that we used both the CDF and E710/E811 high energy experimental 
cross sections at ^/s = 1800 GeV in the ln^(z//m) analysis. Inspection of Fig. |2K) shows that 
at y/s = 1800 GeV, our fit effectively passes below the cross section point of ~ 80 mb (CDF 
collaboration). In particular, to test the sensitivity of our fit to the differences between the 
highest energy accelerator pp cross sections from the Tevatron, we next omitted completely 
the CDF (~ 80 mb) point and refitted the data without it. This fit, also using '^xlmux ~ 6; 
had a renormalized x^/d.f.=1.055, compared to 1.095 with the CDF point included. Since 
you only expect, on average, a Ax^ of ~ 1 for the removal of one point, the removal of the 
CDF point slightly improved the goodness-of-fit. Moreover, the new parameters of the fit 
were only very minimally changed. As an example, the predicted value from the new fit for 
the cross section at a/s = 1800 GeV — without the CDF point — was app = 75.1 ± 0.6 mb, 
where the error is the statistical error due to the errors in the fitted parameters. Conversely, 
the predicted value from Table|2] — which used both the CDF and the E710/E811 point — was 
o'pp = 75.2 ± 0.6 mb, virtually identical. Further, at ^/s = 14 TeV (LHC energy), the fit 
without the CDF point had app = 107.2 ± 1.2, whereas including the CDF point (Table 
= 107.3 ± 1.2. Thus, within errors, there was practically no effect of either 
including or excluding the CDF point. The fit was determined almost exclusively by the 
E710/E811 cross section — presumably because the asymptotic fit was locked into the low 
energy transition energy Uq, thus sampling the rich amount of lower energy data. 

In Table |21 we make high energy predictions of total cross sections and p-values for pp 
and pp scattering — from collider energies up to the high energy regions appropriate to cosmic 
ray air shower experiments. 

We have demonstrated that the duality requirement that high energy cross sections 
smoothly interpolate into the resonance region strongly favors a In s behavior of the asymp- 
totic cross sections for the nucleon-nucleon systems, in agreement with earlier result for 7p 
scattering[6j and up scattering O EI. We conclude that the three hadronic systems, 'jp, Tcp 
and nucleon-nucleon, all have an asymptotic In^ s behavior, thus saturating the Froissart 
bound. 

At 14 TeV, we predict app = 107.3 ±1.1 mb and ppp = 0.132 ±0.001 for the Large Hadron 
Collider — robust predictions that rely critically on the saturation of the Froissart bound. 
Figure 01 shows all available data for both pp and pp, including cosmic ray data previously 
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Table 2: Predictions of high energy pp and pp total cross sections and p- values, from Table ^ 
a ~ W{iy/m^), Ax^^^^ = 6. 



^/s, in GeV 


app, in mb 


Ppp 


(Tpp, in mb 


Ppp 


540 


60.81 ±0.29 


0.137 ±0.002 


60.76 ± 0.29 


0.136 ±0.002 


1,800 


75.19 ±0.55 


0.139 ±0.001 


75.18 ±0.55 


0.139 ±0.001 


14,000 


107.3 ± 1.2 


0.132 ±0.001 


107.3 ± 1.2 


0.132 ±0.001 


50,000 


132.1 ± 1.7 


0.124 ±0.001 


132.1 ± 1.7 


0.124 ±0.001 


100,000 


147.1 ±2.0 


0.120 ±0.001 


147.1 ± 2.0 


0.120 ±0.001 



analyzed by Block, Halzen and Stanev[TT]. It is most striking that the two fitted curves for 
cr„„even, using on the one hand, the ln^(z//m) model of this work and on the other hand, the 
QCD-inspired model of the BHS group^], are virtually indistinguishable over 5 decades of 
cms energy, i.e., in the energy region 3 < ^/s < 10^ GeV. 

4 Part 3: The Glauber analysis for (Jp_air 

The extraction of the pp cross section from the cosmic ray data is a two step process. First, 
one calculates the p-air total cross section from the measured production cross section 



prod cl q—c\ 



p-air "p-air ^p-air ^p-air ■ (-'-3) 

Next, the Glauber method [121 UHl is used to transform the measured value of (y^pl\i^ into a 
proton-proton total cross section a^^; all the necessary steps are calculable in the theory. In 
Eq. (|T3|l the measured cross section for particle production is supplemented with cXpLair and 
CTpZ^lj., the elastic and quasi-elastic cross section, respectively, as calculated by the Glauber 
theory, to obtain the total cross section (Xp-air- The subsequent relation between (Tp'^J^jj. and 

app involves 5, the slope of the forward scattering amplitude for elastic pp scattering, 



where B 



dt \ dt 



, and is shown in Fig. HI which plots B against app, for 5 curves 

t=o 

of different values of a^^^s^i,-- This summarizes the reduction procedure from a^pfl^,. to app [14J. 

The solid curve used the value of B found from the QCD-inspired fit of Block, Halzen 
and Stanev^I], whereas the pp cross section value was found from the fit of Tabled for 
a ~ ln^(i//m) and ^Xima^ ~ 6- The open circle is our value for app at = 72.0 TeV, the 
HiRes energy [in]. 

Unlike the Glauber calculation that was used in references jHj and ^T], we have here 
incorporated inelastic screening into our calculation of eq. (|13p. using a two-channel model 
to approximate diffraction. The screening has the effect of lowering the observed p-air cross 
section, o"p™^i., by about 30 mb at -y/i = 70 GeV. In Fig. El the solid line is a plot of crpL°^i. 
calculated including inelastic screening, against the cms energy y/s. The HiRes point at 72.0 
TeV, the open diamond, is in good agreement with our prediction of CTp^^^^jj.. 

The Fly's Eve|16j and AGASA^T] cosmic ray experiments measure the shower attenu- 
ation length (Am) and not the interaction length of the protons in the atmosphere (Ap_air)- 
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They calculated Ap_air, and thus, crpL°air) from the relation 

A -kX ■ -k^^ (U) 

I\m — ft-^p-air — 1^ prod ' v / 

'^p— air 

where k depends on the rate at which the energy of the primary proton is dissipated into 
electromagnetic shower energy observed in the experiment. The latter effect is parameterized 
in Eq. ()14j) by the parameter k; m is the proton mass and (Tp!°air the inelastic proton-air cross 
section. The values of k in the original publications of ref. P^ ITT] were determined by Monte 
Carlo simulations which did not take into account scaling violations — Fly's Eye used k = 1.6 
and AGASA used k = 1.5. More modern Monte Carlo's such as the SIBYLL simulation [18] 
would give k = 1.2. We have renormalized the k values to fit our curve in Fig. El using 
k = 1.27, in agreement with modern simulations. 

The HiRes point uses a different analysis method [T3] which does not require a k factor, 
being a more absolute measurement of cXp^air- Thus, the agreement of their value of crp™air 
with our prediction is of much more significance — the Fly's Eye and AGASA results are 
only shown because they can easily be made compatible with the predictions of Fig. by 
rescaling to a more modern value, i.e., k = 1.27. 

Finally, we show in Fig. IHl the dependence of the total pp cross section, app on the 
observed p-air production cross section, o"p!°air- 



5 Conclusions 

Using the saturation of the Froissart bound for pp scattering gives us a robust prediction 
of the observed p-air production cross section, crp!°air) when we use the pp extrapolations to 
cosmic ray energies in a Glauber calculation with inelastic screening — thus tying together 
accelerator measurements with precise energy scales to cosmic ray measurements where the 
energy scale has rather large uncertainties. The available cosmic ray measurements of (Tp^fir 
are in reasonable agreement with our predictions. Clearly, more accurate measurements of 
(Tp™air are needed at several different energies, and we urge cosmic ray experimenters to make 
new efforts in these directions. 
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Figure 1: a) A plot of TZ~^, the reciprocal of the factor that multiplies Xmin/'^ found in the fit to 
the sifted data set vs. Axf cut, i.e., ^Ximax- A plot of r^2, the factor whose square multiplies the 

covariant matrix found in the fit to the sifted data set vs. Axf cut, i.e., '^Ximax- These figures are taken 
from ref. Pp. 
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Figure 2: ^^): The fitted total cross sections app and app in mb, vs. ^/s, in GeV;[5l3): The fitted p- values, 
Ppp and ppp vs. -y/s, in GeV, using the 4 constraints of Equations lO, (|10|l . and H12(l . The circles are 
the sieved data for pp scattering and the squares are the sieved data for pp scattering for ^/s > 6 GeV. The 
dash-dotted curve (pp) and the solid curve (jjp) are the fits of Tabled for a ^ ln^{iy/m) and Axf ju^x = ^■ 
The upper sign is for pp and the lower sign is for pp scattering. The short dashed curve (pp) and the long 
dashed curve (pp) are fits of Tabled for a ^ ln(:^/m) and '^Xima.x ~ ^- "^^^ laboratory energy of the 
nucleon is and m is the nucleon mass. 
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Figure 3: The circles are the cross section data for pp scattering and the squares are the cross section data 
for pp scattering, in mb, vs. ^/s, in GeV, for aU of the known accelerator data. The solid curve is the At 
(Tabled a ~ ln^(z//m), ^Xiamx ~ energy data of the crossing-even amplitude, of the form : 

cTnneven = co + ciln (^) + C2ln^ {^j^ + (^)'^ cq and Pv' constrained by eq. © and eq. ifTUIl . 

The dot-dot-dashed curve is the crossing-even amplitude cross section (t„„, from a QCD-inspired fit that fit 
not only the accelerator pp and pp cross sections and p- values, but also fit the AGASA and Fly's Eye cosmic 
ray pp cross sections. 
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Figure 4: B dependence on the total cross section, (Jtotipp)- The five curves are lines of constant cTp^airJ 
of 414, 435, 456, 499 and 542 mb — the central value is the HiRes^Sl value and the others are ±lcr and ±2(t. 
The solid curve is the QCD-inspired^J fit of B against the In^ s fit of atotipp) from Table^for ^XininK ~ ^■ 
The open circle is our value for i/i = 72.0 TeV, the HiRes energy. 
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Figure 5: The p-air production cross section, cTp^^ir' energy, -y/i, in GeV 
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Figure 6: The total pp cross section, app, in mb vs. the p-air production cross section, crpl^^ij., in mb 
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